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Abstract 

The diffusion of a two-dimensional array of particles driven by a constant 
force in the presence of a periodic external potential exhibits a hierarchy of 
dynamical phase transitions when the driving force is varied. This behavior 
can be explained by a simple phenomenological approach which reduces the 
system of strongly interacting particles to weakly interacting quasi-particles 
(kinks). The richness of the strongly coupled system is however not lost 
because, contrary to a single-Brownian particle, the array shows an hysteretic 
behavior even at non-zero temperature. The present investigation can be 
viewed as a first step toward understanding nanotribology. 
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The diffusion of a Brownian particle, driven by an external force and subjected to a 
periodic potential, is a situation which arises in several fields of science such as solid- 
state physics, surface physics 0, chemical physics and even communication theory. It is now 
well understood and provides a simple example of an out-of-equilibrium phase transition, 
between a locked and a running state. The case of many interacting particles is even more 
interesting because collective effects modify qualitatively the picture giving rise, in some 
parameter range, to a dynamical state which is very reminiscent of a traffic jam at the atomic 
scale. We discuss here this behavior in the context of solid state friction because it provides 
a typical example which may be amenable to experimental tests, but the basic ingredients 
required to observe the phenomenon are simple and can be found in many physical systems. 

Understanding the diffusion and mobility of strongly interacting atoms subjected to a 
periodic potential and driven by an external force is a first step toward understanding solid 
friction at the atomic level 0]. For noninteracting atoms the problem would be simple since 
it essentially reduces to the diffusion of a Brownian particle in a periodic potential. Under 
the influence of a dc force F it will preferably diffuse in the direction of the force and in 
average there will be a drift velocity (v) which depends on F. The mobility B is then defined 
a.s B = {v)/F. If, for small forces, B is independent on F (the linear response regime), for 
arbitrary ones a nonlinear response takes place, and the task is to calculate this nonlinear 
mobility. 

The total potential experienced by the Brownian particle is the sum of the periodic 
potential and the potential —Fx due to the driving force, i.e. it corresponds to a corrugated 
plane, with an average slope is determined by F. At small forces the potential has local 
minima, therefore the particle is static and its mobility vanishes. On the contrary for large 
forces there are no stable positions, and the particle slides over the corrugated potential, 
reaching its maximum mobility Bf = (mr^)"^, where m is the mass of the particle and 
rj is the friction coefficient. A simple calculation shows that, denoting by e the height of 
the periodic potential, a the lattice constant, and c a constant depending on the shape of 
potential, the critical force for which the stable positions disappear is F° = ce/a. 
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However, for the underdamped case, the system may also have a running solution, even if 
the minima of the potential do exist. Indeed, because of their momentum the particles may 
overcome the next hill, which is lower than the one from which they are falling due to the 
—Fx contribution to the potential, if the gain in potential energy is greater than the energy 
dissipated during this motion. One finds that this second critical force is = c'rjy/me, 
where c' is a constant which depends on the potential shape. As the particle is either 
locked or running, depending on its initial velocity, the system exhibits bistability and the 
transition between these two states shows hysteresis. However, for a single particle the 
bistability disappears in the presence of an external noise such as thermal fluctuations, no 
matter how small the noise is, because the fluctuations can kick the particle out of the locked 
state by usual thermal activation. Thus, the Brownian motion of a single particle driven by 
an external force shows hysteresis only for zero temperature. 

The case of interacting particles in a periodic potential is a much more difficult problem; 
interesting numerical results were obtained in for high damping, when the time-independent 
Schmoluchowsky equations may be reduced to a one-particle equation with an effective 
on-site potential and then solved numerically by the transfer-integral method [Q. The re- 
sults show that there is a region of highly nonlinear mobility but without a bistability phe- 
nomenon. Besides, recently Persson P] has used molecular dynamics to study a 2D system of 
interacting atoms subjected to a periodic potential. In the underdamped case, he observed a 
dynamical phase transition similar to the one-particle case. Recalling the well-known Aubry 
transition from the pinned state to the freely moving state in the Frenkel-Kontorova model 
with an incommensurate atomic concentration |^, Persson supposed that his results could 
be explained in a similar way. 

In the present Letter, we study the driven motion of interacting atoms in a 2D potential 
for arbitrary damping. First we show that the final state observed by Persson corresponds in 
fact to the sliding motion of atoms in an inclined potential. Then we demonstrate that the 
transition from the locked to the sliding state passes through a hierarchy of hysteretic de- 
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pinning transitions. An important point is that hysteresis persists in the presence of thermal 
fluctuations. Finally, we analytically compute the critical forces and the corresponding mo- 
bilities using a phenomenological approach |^ which treats a system of strongly interacting 
atoms as a system of weakly interacting quasiparticles (kinks). 

Although it is still oversimplified, the generalized Frenkel-Kontorova model that we con- 
sider provides a rather complete description of a layer of atoms adsorbed on a 2D crystalline 
surface when the motions in the vertical dimension are also taken into account. The param- 
eters are chosen for the adsystem Na-W(112) and this model was proposed P] to explain 
the intriguing experimental results obtained for the dependence of the diffusion coefficient 
of strongly interacting adatoms when the concentration varies. 

The model considers the following Langevin equation for the x-coordinate of i-th atom, 



d 

mxi + mrfXi + ^— 



V^{xi) + Vy{yi) + ^mulz^ 



+ J2Voexp{-Po\ri-rj\) = F^^^ + 5Ft\t) , (1) 

and similar equations for the coordinates i/i and Zi. The particles are therefore in a periodic 
rectangular potential, 

Sg (1 + Saf[l - cos{2na/aa)] 
VaW- 2 1 + ^2 _ cos(27r«/a,) ' ^ 

where a is x or y, and parameters s^ = 0.2, Sy = 0.4, = 0.46 eV, Ey = 0.76 eV, = 2.74 A, 
ay = 4.47 A are chosen in relation to the highly anisotropic channeled surface W(112) [Q. 
Atomic mass is m = 1, and we put ujz = 1.84, rj = 0.165 in the corresponding units. The 
exponential interaction law with Vq = 10 eV, f3o = 0.85 corresponds to the repulsion 
of the adatoms at rather high concentration, F^^^ is the external force applied along the 
channels, and 6F is a gaussian random force simulating the interaction with a thermal bath. 
We start with a ground state of the system, then the temperature and later the force are 
increased adiabatically. We then compute the mobility for different values of F and T. 

An important parameter is the atomic concentration 6, corresponding to the ratio be- 
tween the number of particles and the number of available sites M. In this Letter we 



present data at several generic concentrations in order to show various aspects, although the 
general statements are common for the wide range of atomic concentrations. 

Let us first consider the case of = 21/41 modeled by 30 channels, each with = 105 
and M = 205. As this concentration, close to 6' = 1/2, is not a simple 1/q value with q 
integer, the ground state corresponds to large regions where the effective concentration is 
1/2, equidistantly separated by zones of compression p. Since these zones in the standard 
FK model are called kinks (or antikinks in the case of localized expansions), we use this 
terminology here too. The number of kinks is proportional to the difference between 6 and 
the closest value 1/q 0. All the parameters, such as the mass of kink, Peierls-Nabarro 
potential, and energy of creation of kink-antikink pair can be easily evaluated 0. In spite 
of the obvious difference with the mathematical description of soliton behavior in nonlinear 
partial differential equations, we will show that there are surprising similarities and moreover 
that they allow us to explain qualitatively and quantitatively the intriguing behavior of the 
system in a simple way. 

A generic evolution of the mobility as a function of the external force is shown on Fig. 
One first notice a hysteresis and distinguish four different regions when the external force 
is increased. In the very low force range, the mobility is zero. Above a critical force Fj., 
the mobility jumps to a first plateau at Bk- Then there is a second plateau around Bm for 
F > Fpair, and finally for a force higher than Fr, the system reaches the maximum mobility 
Bf. When the force is decreased, the systems jumps directly to the static state for forces 
lower than F^. Let us explain these different states. 

In the first region, the force is too low, and neither the kinks nor the individual atoms 
move. Both are trapped in their wells and the system is in the locked state. When the force 
reaches a critical value F^, the mobility jumps to a nonzero value B^. This state corresponds 
to the kink-running state. Indeed, a careful study of the time dependence of positions of all 
atoms indicates clearly that the compressed zones are moving but not the individual atoms, 
except when a kink passes through their site. A simple analogy with a single particle in the 
periodic potential allows to compute this critical force. Owing to the lattice discreteness, 
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a kink in the FK model moves in the periodic Peierls-Nabarro potential, whose barrier Ek 
could be easily determined 0; in addition, the potential is tilted due to the external force. 
This quasiparticle (kink) will be trapped until the last stable state disappears, i.e. when 
Fk = cEk/a. Besides, as we know the number of kinks 6k in the system, the mobility of the 
kink- running state is Bj. = OkBj. This approach is very successful as shown in Fig. where 
the values and are shown. 

In the kink-running domain, the atoms are static contrary to the kinks since the energy 
barrier e is always greater than the Peierls-Nabarro barrier Ek- Physically this phenomenon 
signifies that it is easier to move a dislocation coherently than to move all the atoms. In 
addition, a detailed study shows a decrease of the relative distance between the kinks. One 
finds a tendency for the kinks to bunch like the cyclists in a "peloton" . The probable reason 
is that, as the kinks are not exact solutions of the system and therefore radiate waves, the 
oscillatory tail of a kink could help the following kink to overcome the Peierls barrier and 
to catch up with the previous one. A second point to notice, is that this motion of kinks 
shows also a hysteresis even for non zero temperature. If one decreases the force when this 
first plateau is reached, the system goes back to the locked state for a critical force lower 
than Ffc. This indicates that kinks have a behavior more complex than a single Brownian 
particle although the plateau at B^ is well predicted by the one-particle picture. 

In the range of force corresponding to the second plateau, at F > Fpair^ additional 
kinks and antikinks could be created since the energy barrier for the nucleation of new pairs 
vanishes. Therefore, we have to take into account not only the geometrical (ground state) 
kinks but also the "force-excited" kinks |^. Fig. ^ presents the time dependence of the 
positions of all atoms. Two different regions can be seen, mobile and immobile ones. The 
finite time interval between the snapshots results in a stroboscopic effect giving a wrong 
impression for atomic trajectories; in order to show one actual trajectory we marked one of 
atoms by black diamonds while others are indicated by unfilled diamonds. The picture is 
very reminiscent of a traffic jam: the particle is trapped into an immobile zone, until being 
first in this region; then the particle ballistically move till the next high-density zone, where 
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it is stopped again. The velocity probability presents a two-bells shape, corresponding to 
static and moving atoms: this is a coexistence regime. 

A careful examination of the B{F) dependence in this range of F shows that the mobility 
is even slightly decreasing with increasing force. Indeed, after the kink-antikink nucleation 
threshold Fpair, the kinks start to bunch into compact groups, as in the context of Josephson 
junction |^. A simple phenomenological theory gives Bm{6) oc Bf {1 — 6)/6, which agrees 
very accurately with the numerical results [|lOl- Moreover, the study shows that the mobile 
subsystem corresponds to a bunch of the antikinks whereas the immobile one corresponds 
to a bunch of the kinks. 

Finally, after the coexistence regime, for high enough force, all atoms are sliding over 
the periodic potential and the mobility reaches its maximum value Bf. When the force is 
reduced, the description of the dynamics in terms of kinks has lost its meaning since there 
is no special organization that subsists. This is why the system goes back directly to the 
locked state at the critical value F^. 

Thus, the system of strongly interacting particles in an anisotropic external potential 
does present a dynamical phase transition when the dc external force is varied. We are able 
to explain the multiple steps by the hierarchy of depinning: first, the geometrical kinks, then 
the "force excited" kinks and finally the atoms. But it is more remarkable that this behavior 
does not disappear with thermal fluctuations as the hysteresis of a single driven particle in 
a periodic potential. Fig. |^ attests that the hysteresis survives for finite temperatures. 

Let us derive an approximate expression for the critical force versus temperature. At 
zero temperature the back transition approximately corresponds for low concentration to 
the external single-particle threshold For nonzero temperatures, the system gets locked 
when the probability for the velocities to be lower than F^/mrj is greater than a threshold Pc- 
As a result, the critical force follows the law 

= + \/^mkBTr]^ erf-^(l - 2P,) = F^ + 6Vt, (3) 

where erf~^ is the inverse of the error function. With F^ ~ 0.144, the solid curve in Fig. |^ 
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shows that this expression scales very accurately with above expression ii S = 0.35. 

Unfortunately, we do not have a complete understanding of evolution of the forward 
transition versus temperature. Because of temperature fluctuations, the particles feel a 
smoother potential and a smaller barrier to overcome; therefore the system jumps to the 
final running state for lower external forces. Numerical results plotted in Fig. |^ for the case 
9 = 21/31 scale with the expression Fr = F,j^ — ^"v/T, if we chose = 0.37 and ^ = 0.44. 
While such a law is valid for other concentrations as well, the parameters significantly depend 
on the concentration contrary to the law for F^. 

The behavior of the ID version of the problem is almost identical to the 2D case. The 
only difference that we have to notice is in the transition to the final sliding state. The exact 
critical force depends slightly on the external conditions, which means that the transitions 
do not occur simultaneously in all the channels. A careful examination of the behavior 
in different channels shows two interesting results. First, when a channel has jumped to 
the sliding state, it enhances the probability for the neighboring channels to jump due to 
interatomic interactions. This first channel corresponds therefore to a nucleation event. 
Indeed, the width of the "river" of neighboring moving channels grows faster than the 
number of independent rivers. The inset of Fig. ^ presents the numerical results which scale 
very well with an exponential law exp [(F - Fq) /AF] with Fq = 0.276 and AF = 0.002. Of 
course, the value of AF attests that it is a very thin effect, but this phenomenon reveals 
an enhanced transition to sliding state due to cooperative effects in the second dimension. 
On the contrary, the transition of a channel to the locked state is almost independent of 
neighboring channels. 

In conclusion, a driven system of interacting particles in a 2-D anisotropic external poten- 
tial exhibits several dynamical transitions succesfully explained by a hierarchy of depinnings. 
This behavior, which is much richer than the behavior of a single Brownian particle illus- 
trates the interesting properties of complex systems. The analysis of the transitions can be 
put in the framework of a single particle theory by introducing collective excitations which 
have particle-like properties although the strong discreteness of the system does not allow 
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us to consider these excitations as solitons. The complex behavior of the multi-particle sys- 
tem is however not lost because, after the transition to the full running state, the collective 
organization is completely destroyed so that the backward transition does not look like the 
upward sequence of transitions. This explains why the complex system can maintain an 
hysteresis at non-zero temperature contrary to the single Brownian particle. This result 
which emphasizes the role of hysteresis at microscopic scale should be related to the main 
role of hysteresis found in solid friction [|lT| . Work along this line is in progress. 
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FIGURES 

FIG. 1. Mobility versus force for a concentration = 21/41 {N = 105, M = 205). The solid 

curve corresponds to increasing force and the dashed curve, to decreasing force. 

FIG. 2. Atomic trajectories in a given channel for the 2D system with 9 = 34/47 at F = 0.24. 
The black diamonds correspond to the trajectory of one atom. 

FIG. 3. The diamonds (triangles) correspond to the position of the transition to the running 
(locked) state for different temperatures in the case = 21/31, while the squares corresponds to 
the transition to the locked state for = 21/41. The threshold was chosen to be S = 0.9Bf. The 
solid and dashed curves correspond to the phenomenological approach discussed in the text. The 
inset shows the width of the "river" (the moving neighboring channels) in the 2D case as a function 
of the driving force and the solid line corresponds to the approximate expression. 
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